Quantum phase shift and neutrino oscillations in a stationary, weak gravitational field 
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We develop a new method for determining the gravitationally induced quantum phase shift for 
a particle propagating in a stationary weak gravitational field, in the framework of the linear ap- 
proximation. This method is based on the properties of the Synge's world-function. The result is 
applied to the neutrino oscillations in the field of an isolated, axisymmetric rotating body. The 
■ gravitational field is described using the linearized part of the post-Newtonian metric yielded by 

f***) ' the standard Nordtvedt-Will formalism. An explicit calculation of the phase shift is performed for 

, neutrinos produced inside a homogeneous nearly spherical rotating body. Our general formulae are 

£\j . valid even for nonradial propagations. The gravitational effects found here involve only the PPN 



S3 
in 



o 



parameter 7. 



I. INTRODUCTION 



In the last fewyears, a series of papers has been devoted to the effects of gravitation on the neutrino oscillation phases 
\Q ' (see Refs. jj], ||, 0, ^] and references therein). As far as we know, all the proposed derivations of the neutrino oscillation 
formula within the WKB approximation use the explicit determination of the timelike geodesies corresponding to 
freely falling particles. For this reason, even in the simple case of the Schwarzschild metric treated in the linear 
approximation, the calculations become heavy for nonradial propagations and the results remain limited to the case 
where the neutrinos are created outside the matter generating the gravitational field. Our purpose in the present paper 
is to provide a method based on the world- function developed by Synge [^| , which does not require the determination 
of the geodesies. This method, which has recently proved its worth for electromagnetic signals ||, can be easily 
applied to neutrinos and more generally to all kinds of particles, as long as we may ignore the spin contributions [Q. 



We begin by giving a derivation of the quantum mechanical phase valid for any particle freely propagating in a 
■ given stationary, weak gravitational field. We determine both the travel time and the accumulated phase between 
rr| two points of space-time in the limit of the semi-classical approximation, where the particle is assumed to follow its 
!>" ■ classical path. Next we apply the results to the neutrino oscillation phases. We adopt the general assumptions of 
KJ{ ' Bhattacharya et al. || concerning the interference of two mass eigenstates observed at the same point of space-time. 
For applications, we consider the linear terms in G of the post-Newtonian metric given by the standard Nordtvedt-Will 
formalism |^|. This kind of approximation enables us to carry out the calculations for neutrinos travelling across the 
gravitational field of an isolated, axisymmetric rotating body. We obtain explicit formulae for a homogeneous, nearly 
spherical rotating body, including the case of nonradial propagations. 

In this paper, G is the Newtonian gravitational constant and c is the speed of light in a vacuum. We suppose that 
space-time is covered by a global quasi Cartesian coordinates system (x M ). We put x° — ct, t being a time coordinate 
and we use the notation x = {x l ). The signature adopted for the metric tensor g^ v is (H ). 
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II. PHASE SHIFT IN A WEAK STATIONARY GRAVITATIONAL FIELD 
A. Phase shift and world-function 

Let us consider a free particle with proper mass m propagating from a point x A to a point x B . The quantum 
mechanical phase of this particule accumulated between x a and xb is given by 

Trie 

$AB = -T-SAB , (1) 

a 

where sab is the path length of the timelike geodesic Tab joining x a and xb ||- Denoting by A the unique affine 
parameter along Tab which satisfies the boundary conditions Xa = and Xb = 1, we can write 

s AB =2Q{x A ,x B ), (2) 

where Q(x a ,xb) is the so-called world-function, defined as 

1 r 1 dr^ dr u 

Sl(x A ,x B ) = -J o 9tny (x a (X)) ——dX, (3) 

the integral being taken over Tab- Consequently, the accumulated phase <&ab may be written as 



= —y/2n(x A ,X B ). (4) 

It follows from Eq. (Q) that the knowledge of the world-function corresponding to a given metric enables to 
determine the accumulated phase $ab- In addition, the world- function enables to perform the calculation of the 
energy-momentum vector p a = mcu a of the particle at xa and xb (u a denotes the unit tangent vector dx a jds along 
Tab)- Indeed, it is shown in Ref. Q that the covariant components of u at xa and xb are given by the relations 

{u a ) A = 7= , {U a )B = —7 =, , TT-^ ■ 5 

y/2Q(x A ,x B ) ® x a ^/2n(x A ,x B ) 9x% 

These relations are very useful, especially in the case where the metric does not depend on one or several coordinates, 
as we shall see in the next subsection. 

B. Phase shift in a stationary space-time 

Throughout this paper, we restrict our attention to stationary gravitational fields and we choose x a so that does 
not depend on x°. The world- function has then the form fl(x° B — x° A , x A , xb) and the component uq of the tangent 
vector to the classical trajectory is a constant of the motion, i.e. 

i*o = E/mc 2 , E = const. (6) 

When the particle is passing through a point x, the local energy as measured by a stationary observer staying at x 
is Ei oc (x) = E/ ^/go (x), so the constant E may be considered as the asymptotic energy of the particle at spatial 
infinity. It follows from Eqs. (||) that 



dfl(x° B - x° A , x A , xb) _ E 
^2n{x B ~x A ,x~ A ~^ B ) d(x B ~ x° A ) 



u - rr^r n n =^ HTZn Zol ~ 33 ■ i 7 ) 



Eliminating 2Q(xa,%b) between Eq. (Q) and Eq. ([?]), we obtain 

m 2 c 3 dn{x° B -x° A ,x A ,x B ) 
* AB ~^E d(x B -x A ) ■ (8) 

Then, eliminating x B — x A between Eqs. ^ and (|^), it will be possible - at least in principle - to carry out the 
calculation of the accumulated phase <&ab as a function of x A , x B and E. 
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C. Linear weak field approximation 

For subsequent applications, we suppose that the gravitational field is weak. So, the metric may be written in the 
form 



9V = Vim + !l 



(9) 

where r)^ = diag(l, —1,-1,-1) and where the quantities are considered as perturbation terms of the first order 
in the Newtonian constant G. The word-function Q(xa, x b ) may be expanded in powers of the gravitational constant 



n(x A ,x B ) = n {0) (xA,x B ) + n {l \x A ,x B ) + o(h 2 ) , 

where CI^°'(xai x b) is the world-function in Minkowski space-time 

^\ XA , XB ) = \[{ x %~ XA y-R AB \ , 

with Ra B = \x b — xa\, and Q^(xa,x b ) is given by 

Q,W(xa,x b ) = \{x B -x A )(x B -x A ) [ h^ u (x( ) (A))dA, 



(10) 
(11) 

(12) 

the integration being performed along the path of parametrized equations x^(X) = [x B — xa)X+xa- Since space-time 
is assumed to be stationary, Eqs. (|l0|), ( pi] ) and (1^) yield 



Q.{x a ,x b ) = -{x b ~ Xa Y[1 + 2A{x a ,x b )] + { Xb -x a )R ab B{x Ai x b ) 
-Was [l-2C(x A ,x B )}+0(h 2 ), 



(13) 



in which we have set 



A(x A ,x B ) = -j hoo(x( )(X))dX, 
B(x A ,x B )= ^ x% b- x a) f ho .( X{0) (\)) dX . 



C(X A ;X B ) 



Ra B Jo 
( X B — x A )(x B — x A ) f 

tR\ B 



(14) 



hij{x(o){X))dX. 



with 



a; (o)(A) = (x B - x A )X + x A 



(15) 



For the sake of brevity, we shall henceforth omit the symbol 0(h 2 ) standing for the post-post-Minkowskian terms. 
Differentiating the r.h.s. of Eq. (|l3| ) with respect to x B — X ° A and taking Eq. (||) into account, we obtain for 



®ai 



m 2 c 3 
HE 



[(x° B - .<) (1 + 2A{x A , x B )) + RabB(x a ,x b )] 



(16) 



Now, substituting for f2 from its expression ( ^3|) into Eq. (|]) and solving the equation thus obtained for X ° B — x Al we 
find 



x B — x A — R A 




2 A 

m c 



1 ~ I 1 ~ E 2 _ m 2 c 4 ) A i x A,x B ) - C{x Al x B ) 



- B(x A ,x B ) 



} • (17) 



Finally, inserting into (|l6|), we get an expression for <&a B as follows 



§A B = 



Rae m c 



2,,:s 



Ti \JE 2 - m 2 c 



1 



E 2 



E 2 



-A(x A ,x B ) - C(x A ,x B ] 



(18) 



The formulae (|17|) and ( |18[ ) giving respectively the time delay and the accumulated phase between xa and Xb are 
valid whatever the speed of the particle. 
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III. NEUTRINO OSCILLATIONS IN THE FIELD OF AN ISOLATED, AXISYMMETRIC ROTATING 

BODY 

A. General formulae 



Let us apply the previous results to the neutrino oscillation phases. We suppose that the neutrinos are produced 
at fixed local energy Ei oc in a weak flavour eigenstate which is a linear superposition of two mass eigenstates mi and 
77i2, with 777,2 < TTii. Moreover, we suppose that the source of neutrinos is located at a fixed point xa and that the 
interference is observed at the same final space-time point (x%, Xb)- As a consequence, the eigenstates of masses mi 
and 777,2 are respectively emitted at the times (x A )\ and (x A )2, with (x° A )i < {x° A )2- It is generally considered Q that 
the phase difference A<I> = $2 — $1 between the two mass eigenstates at the point xb is given by 



A*=|[M) 2 



(x° A )i] + (<Pab)2 - (®ab)i 



which can be rewritten in the form 



E 
he 



X% - (x A )t (x% - (4) 2 )] + ($ AB ) 2 - ($ABh 



Using (pj 



A$ = - 



and (§1 

Rab 



it is easily seen that the phase difference A$ is given by 



(m\ — m^c 3 



n y/E 2 - m{c A + y/W 



1 



E 2 



,2,4 



VE 2 



,2,4 



A(x A ,x B ) - C(x A ,x B ) 



(19) 



(20) 



(21) 



We note that the quantity B(xa,xb) does not appear in expression (|T|). According to (|i"4|), this means that the 
components goi of the metric do not contribute to the phase difference A$ in the limit of the weak- field, linearized 
approximation. 

We can apply the formula (|2l| ) when the stationary field is treated in the slow-motion, post-Newtonian limit of 
metric theories of gravity, provided only the metric perturbations of order 1/c 2 and 1/c 3 are taken into account. We 
use the Nordtvedt-Will PPN formalism and we suppose that the coordinates (x M ) are a standard post-Newtonian 
gauge S. Since the potentials hoi do not influence the phase oscillations, the only relevant potentials are then 



^00 = — nU 



(22) 



• 0(4) , h l3 = — Mj + 0(4) , 

c~ C 

where U is the gravitational potential. From Eqs. (|l4|), we deduce 

A(x A ,x B ) = J f7(cc(o)(A))c2A, C(x A , x B ) = J U{x {0) (X))d\ , 

where a;( )(A) is given by (|l5|). Substituting for A(xa,xb) and C(xa,xb) from Eqs. ( p3| ) into Eq. (|2l|), we obtain 
the general formula for the phase difference 



(23) 



A$ = - 



Rab 



(m 2 — m^c 3 



^JE 2 - m 2 c A + ^JE 2 




E 2 



7 



^JE 2 - m 2 c 4 ^JE 2 



U(x (Q) (X))d\ . (24) 



This formula may be applied to neutrinos propagating in the field of an isolated, axisymmetric rotating body. It 
is immediately seen that the angular momentum of the rotating body has no influence on the neutrino oscillations in 
the weak-field limit and that there exists no coupling with the velocity v r of the center of mass of the body relative 
to the universe rest frame. This last feature is natural since 7 is the only PPN parameter involved in the expression 
ofA$. 

Specializing to the case where mic 2 <ti E and 7772c 2 <C E and expanding the formula (|24|) in series of powers of 
mic 2 / E and of 7772c 2 / 'E, we find 



A$ = - 



Rab (ml 



m 2 )c 3 



2/i 



E 



1 



7- 



1 



U(x, 0) (\))d\ 



( m i 



4E 2 



1 



7 



U(x/ )(X))dX 



This development shows that in general relativity, there is no contribution of the Newtonian potential U up to the 
order m 2 — m\ since 7 = 1. This conclusion holds equally well for a propagation inside the body as for a propagation 
outside the body. Thus we significantly extend a result previously obtained for the exterior Schwarzschild metric 
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B. Neutrinos in the field of a nearly spherically symmetric rotating body 



According to and (25), an explicit calculation of the phase difference A$ requires only to perform the integration 
of the Newtonian potential U along the line defined by (jig). As an example, we may consider the gravitational field 
generated by an isolated, nearly spherically symmetric rotating body. For the sake of brevity, we neglect the multipole 
structure of the gravitational field, which means that the body is characterized by its mass M and by an external 
radius r . We shall examine two cases. 

The geodesic path joining xa and xb is entirely outside the body. In this case, we have U(r) — GM/r, which yields 

1 TTt t\\\j\ GM r A +r B + R AB , 

U(x (0) (X ))d\ = - — In — , 26 

Rab r A +r B - Rab 

where ta = | xa I an d r B =\xb \ ■ 

Substituting ( gtj ) into ( p5j ) yields the explicit expression of the phase difference for nonradial geodesies outside the 
body. For radial geodesies, the formula ( pq ) reduces to 

U(x (0) ))d\= GM In — . (27) 

r B - TA TA 

Inserting (^?j) into ( f25j ) and putting 7 = 1, we recover up to the order m\ — m\ a result previously given in Ref. |l(J 
for the radial geodesies in the Schwarzschild metric (this effect was also calculated by the authors of Ref. ||) , but the 
terms of order ra\ — m\ in their Eq. (13) must be divided by 2). More generally, substituting (1 + u>)/(2 + uj) for 7 
into (pjj), we get the expression of A$ within the frame of the Bergmann- Wagoner tensor-scalar theories of gravity 
(see, e.g., || and references therein). Thus we extend a result previously found (T^J at the order m\ — m\ for radial 
geodesies. 

The point of emission xa is located inside the massive body and the point of reception x B is located outside. In 
order to simplify the evaluation of the integrals in the r.h.s. of (p5|), it is useful to introduce the new affine parameter 
A defined by 

A = A+ r * ^ 2 Rab , (28) 



1Rab 



the range of A being A^ < A < As, with 



r 2 _ 2 _ r,2 2 _ 2 , n>2 

T r B r A n AB T r B A' U AB ( 0Q \ 

Ltx AB zrL AB 

It is easily seen that the value Ao corresponding to the intersection of the segment connecting xa and x B with the 
sphere r — ro is given by 



Rai 



(30) 



where r p is the impact parameter determined by 



l4 (XA ' XB? -\(r\ + r%) - \r\ b - {r \ ~f A? . (31) 



We have 



d2 — o \' A 1 ' B A ''AB at>2 

R AB 2 4 AR AB 



1 /.An /-A 



U(x [0) (X))dX= l_ U l (x [0) (X))dX+ I GM T , dX, (32) 

J\ A J\o \ x (0){X)\ 

where Ui is the potential inside the body. This relation holds for any mass density law compatible with the spherical 
symmetry. For the sake of simplicity, let us assume that the body is homogeneous. Then, Ui is given by 

GM 



Ui(r) = -— 3--y) , 0<r<r o , (33) 

Zrt) \ Try 
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and a straightforward calculation yields 

GM 

U(x {0) (\))d\ = — 



r A +r„ + e^/rfi-r^/ri -r] 
6^ 



In 



rB 



Rab 



2R 



AB 



\ 



rl-rl 



J 



(34) 



where e = 1 if R AB < r 2 B - 
containing e as a factor vanish). 



and e = -1 if R 2 AB > r B - r 2 (if R 2 AB = 



b 1 p 



r 2 ,, then r A — r 2 — and all terms 



IV. CONCLUSION 



It is clear that the world-function constitutes a very powerful tool for finding the quantum phase shift for particles 
freely propagating in a stationary weak gravitational field. The method works equally well for propagations inside 
matter as for propagations outside matter, provided the spin effects may be neglected. 

Using this method, we have obtained a set of formulae giving the gravitationally induced neutrino oscillations in the 
field of an isolated, axisymmetric rotating body. The Newtonian potential is the only gravitational quantity occuring 
in these formulae, which implies that the neutrino oscillations are not influenced by the internal angular momentum 
of the source of the gravitational field. 

Owing to the simplicity of the method, it has been possible to perform explicit calculations for neutrinos propagating 
in the field of a homogeneous, nearly spherically symmetric rotating body. Let us emphasize that these calculations 
are straightforward even for neutrinos emitted inside the body and travelling along nonradial geodesies. 
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